Order Parameters and Phase Diagram of Multiferroic RMn 2 Os 
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The generic magnetic phase diagram of multiferroic RM^Os (with R=Y, Ho, Tb, Er, Tm), 
which allows different sequences of ordered magnetic structures for different R's and different control 
parameters, is described using order parameters which explicitly incorporate the magnetic symmetry. 
A phenomenological magneto-electric coupling is used to explain why some of these magnetic phases 
are also ferroelectric. Several new experiments, which can test this theory, are proposed. 
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There has been much recent interest in multiferroics, 
which display simultaneous magnetic and ferroelectric 
(FE) ordering US HQ. In particular, the orthorhombic 
family RM112O5 (RMO), where R is a rare earth, exhibits 
interesting sequences of magnetic density wave orderings, 
with varyingwave vector q, and some of these phases 
are also FE @, @, 0, H, Q ■ In all these phases one has 
q y = 0, while \q x - || < 0.02 and \q z - || < 0.02. Cool- 
ing from the paramagnetic (PM) phase, one first enters a 
phase in which both q x and q z arc incommensurate. We 
call this phase Hi (1= "incommensurate" , and the sub- 
script will be explained below; some experimental papers 
call this phase 2DIC). For R=Y , Er @| and Tm [| , fur- 
ther cooling yields transitions into a phase which we call 
IC2 (also called 1DIC), where q x is still incommensurate, 
while q z = \ (C= "commensurate"), then into a "CC" 
phase (also called CM), with q = (|, 0, j), and finally 
into a phase where both q x and q z are incommensurate 
again ("II 2 ", or LTI-2DIC). R=Ho Q and Tb go di- 
rectly from Hi to CC. For R=Er, the low temperature 
(T) phase has q x = h, while q z is incommensurate ( "CI" , 
or LTI-1DIC). While the phases IC 2 and CC exhibit a 
FE moment P along the y (b) axis, such a moment ap- 



prcdictions, which can be tested experimentally. 



pears in only some of the observed low T phases (lfj, [11 1 . 
Up to now, the microscopic theories of these systems are 
controversial, and a phenomenological description which 
provides a unified explanation of this complicated phase 
behavior does not exist. The present paper rectifies this 
situation, and provides a basis for analyses of other mul- 
tiferroics with large unit cells. 

Although group theory has been applied to neutron 
diffraction data from magnetic materials [l2^ . its impli- 
cations for multiferroics have not been fully exploited 
until the definitive analyses of NI3V2O8 and TbMnO.3 
@, 0, O , HI • Following the same approach, we iden- 
tify the order parameters (OP's) allowed by symmetry 
[la ] and find the generic phase diagram for RMO sys- 
tems (Fig. [T]), which allows for the observed sequences 
of phases. The theory also explains (a) which phases are 
simultaneously magnetic and ferroelectric, (b) the occur- 
rence of two distinct spin structures in neutron diffraction 
studies of the CC phase 0, 0], and makes several new 
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FIG. 1: LHS (supplied separately): Schematic 3D phase dia- 
gram for q x (q z ) near i (|). J x and J z are parameters which 
control q x and q z . The red surface separates PM and Hi . Be- 
low the blue surface one has q z — \, in phases ICi and IC2. 
The green surfaces represent Hi^Il2 and ICi— >IC2. Below 
the orange surfaces q x — ^, in phases CI or CC. RHS: a cut 
at constant q x 7^ 1/2. The dashed and dotted lines represent 
proposed trajectories for specific RMO's, as T is varied. 



The PM unit cell of the RMO's contains 4 Mn 3+ , 4 
Mn 4+ and 4 R 3+ ions. Denoting these ions by r(= 
1, . . . , 12), and the corresponding Fourier transforms of 
the a-spin-components by S a (q, t) , the quadratic terms 
in the Landau free energy Fm are 



Fm,2 = 



- E 

2 v 

q,a,p;r,r 
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In principle one would diagonalize the (36 x 36) in- 
verse susceptibility matrix Xatt( c k T ' T ') (determined by 
the various magnetic interactions). As T is lowered, the 
first phase to order corresponds to the eigenvalue which 
approaches zero first. The degeneracy of this eigenvalue 
has two origins: first, all the n q wave vectors in the star of 
symmetry-related optimalwave vectors q's have the same 
eigenvalue. Second, each of these q's is associated with 
irreducible representations (irrep's) T of the PM symme- 
try group of q (the 'little group') Excepting acci- 
dental degeneracy, a continuous transition from the PM 
phase involves only a single irrep. If this critical irrep is 
d dimensional (dD), then this eigenvalue is dn 9 -fold de- 
generate and this manifold is described by dn q real OP's, 
or driq/2 complex ones. Each complex OP represents the 
amplitude and the phase of the spin ordering cigenfunc- 
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tion, {^(q, r)}. The symmetry of the cigcnfunction is 
associated with the irrep and is inherited by the OP's. 

For each RMO, the optimal wave vector q is deter- 
mined by its specific material (e. g. the exchange and 
anisotropy energies) and experimental (e. g. pressure, 
magnetic field) parameters. We represent these control 
parameters by their combinations, denoted J x and J z , 
which fix the values of q x and q z , respectively. Figure [1] 
shows the phase diagram of the RMO's in terms of J x 
and J z . Following experiments, we fix q y = 0. We start 
with the case q x ^ \ (with q z near 4). For each such 
q, the 'little' group contains only unity and m y , which 
maps [x, y, z) into (x+^, y+^, z). This group has two ID 
irreps, r a and Tf,, with complex OP's er a (q) and cr fc (q). 
Inversion symmetry 2 then implies non-trivial relations 
between the 5 a (q)'s and the Sp(— q)'s (which have the 
same eigenvalue), reducing the number of independent 
parameters. This should ease future accurate analyses of 
the neutron data. For r a , symmetry implies 



m„cr 



(2) 



and similarly for T b (p depends on the origin) [20| . 

For q x ^ i, the star of q contains four wave vec- 
tors, namely, q± = [q x , 0, ±q z ) and — q±. Therefore, 

we introduce two complex OP's, <r+ = a a (q+) and 
a = "aCl- )j associated with irrep T a and similarly 
for r 6 . Then, F M>2 = E s = a ,t( T ~ T C,s)[Wf? + 
Rejecting accidental degeneracy, we set Tc, a > Tc,b and 
identify the 2DIC phase with our Hi phase, associated 
with a single irrep (the subscript 1 refers to the number 
of irreps), represented by the u^'s. The transition PM 
— > Hi occurs at T = Tc\ a , represented by the top (red) 
surface in Fig. [TJ Which OP's actually order depends on 
the quartic terms in the free energy. For q x ^ these 
include 
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+S G (W aa [at (a~y] 2 + cc) S [G - (0, 0, 4q z )] , (3) 

where G is a reciprocal lattice vector. For q z ^ \ and 
T < T G , a one has \<r+\ = \a~\ > if U a < 0, and only 
one of the OP's orders otherwise. For q z near j, the 
Umklapp term with W aa locks q z to j, in a phase called 
ICi. Within Landau theory, this happens below a first 
order surface (blue in Fig. [1]), parabolic in J z . 

As T is reduced, more quartic terms need to 
be considered, notably W J2 m =±{[ (7 a( c im ) )<7 b {q L m ) )*? + 
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Assuming that and q^ ; are al 



most the same, this term locks the optimal q^ to q±', 
at some T slightly below Tc, a , where a b has not yet or- 
dered (without involving a phase transition.) Accord- 
ingly, we no longer keep the superscripts (a, b) on the 
q's. As T is further reduced, the tendency of the spins 
to have fixed length (rather than oscillate sinusoidally) 



[H, 21 1 may cause a second continuous transition, into the 
phase II2 (or IC2); where both a a and a b are nonzero. As 
shown below, this transition (green surface in Fig. [I} oc- 
curs at a temperature which is parabolic in J x . 

We next discuss the special case q x = \ ( or Jx = 
J x , c , at the back of the 3D diagram in Fig. [TJ. For 
q = (i, 0, g z ), the little symmetry group changes: it now 
contains the additional glide operation m x [which maps 
(x, y, z) into (— x + + i, z)]. This group has only one 
2D irrep [l4, 22 1, with two degenerate complex OP's, a 1 
and a 2 , and corresponding eigenvectors as listed in Table 
XVI of Ref. E3 Hi]. These OP's transform as [Hi 



where ( n = (-1) 



Cn&3- n , I&r, 
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= 1,2. Cooling from the PM 
phase, exactly at J x = J x>c and q z ^ j, one first goes 
into the CI phase, with the free energy 

F M = {T-T c )[\a^ + \a 2 \ 2 ]+u[\a^ + \a 2 \ 2 f 
+ W c \ 

On further cooling, additional Umklapp terms cause a 
first order transition into the CC phase where q = 
(^,0, j). This lock-in happens under a parabola, which 
connects to the blue parabolas which appear for q x ^ \- 
We next vary J x away from J x . c . The RHS of Fig. 
UJ shows a cut of the 3D phase diagram, at fixed J x — 
J XtC + AJ. For small A J, mirror symmetry, q x — > — q x , 
implies that the inverse susceptibility has two branches 
of eigenvalues given by xjjj 1 (&) — T — Tc + ak' x ± bk x AJ, 
where k x = 1/2 — q x and a and b are constants. At 
quadratic order (i. e. using Fm,2)i this implies that 
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T c ,b = 2bk x 0) AJ, where k ( x ] = bAJ/(2a) mini- 
mizes Xi^ilx)- This gives rise to the (green) Hi — >Il2 
phase boundary, Tc. a — Tc.b oc (AJ) 2 . 

For q x close to |, further cooling may lock it to |, 
due to Umklapp terms. This may generate transitions 
into the CI or the CC phase, for T's below the orange 
surfaces in Fig. [TJ A detailed analysis shows that these 
surfaces are also parabolic in A J. The actual sequence of 
transitions then depends on which parabola is narrower. 
In Fig. [ljwe show the case when the green parabolas are 
broader than the orange ones. In this case, the orange 
surface represents II2— >CI and IC2^CC. In the opposite 
case, the phases II2 and IC2 never appear. As shown 
in Fig. [TJ both parabolas are shifted upwards below the 
blue surface, where q z — j, due to Umklapp terms. 

Equations (|T][5j lead to a natural interpretation of neu- 
tron scattering results for the CC phase in YMO. Figure 
[5] shows the Mn 3+ a-b plane spin components in the CC 
phase of YMO, from the data of Refs. d andGj. 
These two structures are obviously similar, and one might 
ask what symmetry (if any) relates them [2II ] . Since the 
structure on the left (right) is even (odd) under the glide 
operation m x , we conclude that the structure on the left 
(right) has er 2 = {<7 l = 0). Going between these two 
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structures corresponds to a rotation in OP space; the in- 
plane spin components belong to distinct but equivalent 
structures. Since either a 1 = or a 2 = 0, we conclude 
that in Eq. flS|), the net coeficient of \cr 1 a 2 \ 2 (Wc — 4|Vc| 
plus the additional Umklapp terms) is positive, prevent- 
ing both OP's from ordering simultaneously [14[. 
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FIG. 2: (Color online) Schematic diagram of the a and b 
components of the Mn 3+ spins in a single a-b plane of YMO 
for the CC phase. The glide m x consists of a mirror plane M 
at x = a/4 followed by a translation b/2 along y. LHS: the 
structure given in Table III of Ref. (with the c-components 
not shown). RHS: the structure given in Fig. 2 [24| of Ref. 
[l8l (who reported zero c-components of spin) . 

For the phenomenological description of the multifer- 
roicity, the total free energy is Fme = -Fm + ^P 2 ^ 1 +Vint, 
where e is the dielectric susceptibility and VS n t the magne- 
toelectric (ME) interaction. Although such interactions 
can also generate a spatially non- uniform P, here wc dis- 
cuss only the uniform case. We again start with the II 
and IC phases, where q x ^ h. To lowest-order, wave vec- 
tor conservation and time-reversal invariance give [H, 27 1 



(0) 



s,t=a,6q=q± 7 



The terms with s = t vanish because they arc odd under 
X. For the Hi (ICi) phase, only a a is nonzero, and there- 
fore P = 0. To have P / with q x ^ | we must have 
the superposition of two irreps, and this happens only in 
the II2 or the IC2 phases. In these phases, we have 



hit 



E E[ f 

q=q± 7 



yo- a (q) cr b(-q) + cc]p 7 



(7) 



and invariance under I requires that r 7 is real. From Eq. 

<J a a b i s °dd under m y . For V5nt to be invariant under 
m yi P 1 must be odd under m y : symmetry forces P to be 
along y(b), as observed ( Hig her order ME interactions 
weakly violate this result [20|]). . 

In the CC phase, Eq. ([6]) is invariant under the sym- 
metry operations of Eq. ([4]) only if [lj| 



const.x[|(T 1 (q)| 2 -|t7 2 (q)| a ]P v 



(8) 



Note that Eqs. ([7]) and (JSJ) apply whether the micro- 
scopic ME interactions are due to exchange striction [28[ 
or to charge orderin g |2 9l|. Thus, P must lie along y(b) 
also in the CC phase [30(. Within mean field theory, Pb is 
proportional to |(er 2 )|, as is the intensity of the magnetic 
Bragg peaks. This is confirmed 17j in RbFe(Mo04)2 
[which is also described by Eq. ©] and also apparently 



for ErMO by Ref. [HI (32|. Since the CC phase is fer- 
roelectric, the fourth order terms in Eq. (0 (plus the 
Umklapp terms) must select <J x a 2 = 0, which we deduced 
from Fig. [21 (The alternative would imply \a x \ = |<7 2 |, 
hence P = 0.) In fact, the selection of which OP, a x or 
cr 2 , is nonzero is a result of broken symmetry. An electric 
field along y(b) would order P y , and then Eq. ([8]) would 
select either a 1 or a 2 , depending on the sign of the field. 
Therefore we suggest that the sample should be cooled 
into the FE phase in the presence of a small electric field 
along y. Depending on the sign of the electric field one 
should get either the left-hand or the right-hand panel of 
Fig. 2 [Si. 

Equation has further implications. First, near the 
P— >IIi transition, a leading fluctuation expansion yields 
Ae oc (P b 2 ) cx \{o~a\){o~1)\- Since only a a becomes critical 
there, we expect singularities in e which behave as the en- 
ergy (IT-Tdl 1 "") and (for T < T C1 ) as the square of the 
OP ((T C1 -TW, but with n = 4 exponents [3J]. Indeed, 
experiments [35| show a break in slope at Tq\ , apparently 
confirming this prediction. In addition, this anomaly in 
the zero frequency dielectric function e{uo = 0) reflects the 
emergence of a resonance in e(w), due to electromagnons 
[36| . Second, in the Hi phase (a a ) ^ 0, so that Vi n t be- 
comes — 2rb<s[{o' a )al\Py. This bilinear coupling between 
P y and 3[(<7 a )0ft] has several implications on the criti- 
cal behavior near the Hi— >II 2 
discovered in some new RMO 



should this transition be 
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Finally, we associate the different RMO's with trajec- 
tories on our phase diagram. Since ErMO Q , TmMO 0] , 
and YMO [B| exhibit ferroelectricity in the phase denoted 
1DIC, wc must identify this phase with our IC 2 phase, 
where both a a and a b order. For these materials, the 
experimental path in parameter space apparently goes 
from Hi via IC 2 into the CC phase, as indicated by the 
dashed lines in the RHS of Fig. [1] (these lines have small 
slopes, since the experimental optimal q varies with T: 
J x and J z depend on T due to other degrees of freedom). 
The last term in Eq. (J3j> implies that q z locks to j only 
if both q + and q_ appear in the ICi and IC 2 phases. If 
2|Woo| > U a > 0, then the two q's first appear as the ICi 
phase is entered. If U a < 0, then both wave vectors would 
have already condensed simultaneously in the Hi phase. 
It would be interesting to determine which scenario ac- 
tually occurs. Since the ME interaction is significant, we 
suggest to apply an electric field parallel to one of the 
q's, and check whether in the Hi phase the two q's arise 
in separate domains or coexist within a single domain, 
following the logic of Ref. [38|. In contrast to the above 
RMO's, HoMO @] or TbMO go directly from 2DIC 
to CM, as along the dotted line in the RHS of Fig. Q] 
Both sequences are thus allowed by our theory. 

The Landau theory is probably less useful at lower T: 
the low T phases depend on the details of the magnetic 
interactions, and higher order terms in Fm should be 
included. Such terms could turn the (orange or blue) 
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surface bounding the CC phase backwards, thus allow- 
ing transitions back into the paraelectric Hi phase, the 
weakly FE phases II 2 or IC2 or the FE phase CI. Also, the 
trajectory describing each material need not be straight 
(thick dashed line in Fig. [T|). A parabolic line, like the 
thin dashed line, would yield a transition from CC to II2 
(or even to Hi) with decreasing T. In fact, in ErMO Q 
the LTI phase seems to have q x = A , which identifies this 
phase with our CI phase. Thus, the observed LTI phase 
could be any of the phases on the other side of the CC re- 
gion, paraelectric or weakly ferroelectric. The effects of a 
magnetic field can be explained as follows: the field gen- 
erates magnetic moments on the R ion (even above their 
ordering temperature). Since these ions couple to the 
Mn ions, their moment results in changes in the effective 
Mn-Mn interactions, thus changing the 'control param- 
eters' and the optimal q. Apparently, this often moves 
the material towards the CM regime, resulting in a tran- 
sition from the low T phase (Hi or II2) back into the CC 
phase . Similar effects happen due to pressure [35| . 

Neutron diffraction measurements in a magnetic field and 
pressure could help resolve these scenarios. 

In summary: we have developed a phase diagram to 
explain the multifcrroic behavior of the family of RMO 
systems and have proposed several experiments to ex- 
plore the unusual symmetries of these systems. 
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